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Abstract — The calculation of network reliability in a proba- 
bilistic context has long been an issue of practical and academic 
importance. Conventional approaches (determination of bounds, 
sums of disjoint products algorithms, Monte Carlo evaluations, 
studies of the reliability polynomials, etc.) only provide approxi- 
mations when the network's size increases, even when nodes do 
not fail and all edges have the same reliability p. 

We consider here a directed, generic graph of arbitrary size 
mimicking real-life long-haul communication networks, and give 
the exact, analytical solution for the two-terminal reliability. 
This solution involves a product of transfer matrices, in which 
individual reliabilities of edges and nodes are taken into account. 
The special case of identical edge and node reliabilities (p and 
p, respectively) is addressed. We consider a case study based 
on a commonly-used configuration, and assess the influence 
of the edges being directed (or not) on various measures of 
network performance. While the two-terminal reliability, the 
failure frequency and the failure rate of the connection are 
quite similar, the locations of complex zeros of the two-terminal 
reliability polynomials exhibit strong differences, and various 
structure transitions at specific values of p. 

The present work could be extended to provide a catalog of 
exactly solvable networks in terms of reliability, which could be 
useful as building blocks for new and improved bounds, as well 
as benchmarks, in the general case. 

I. Introduction 

Network reliability has long been a practical issue, and 
will remain so for years, since networks have entered an 
era of Quahty of Service (QoS). IP networks, mobile phone 
networks, transportation networks, electrical power networks, 
etc., have become "commodities." Connection availability rates 
of 99.999% are a goal for telecommunication network opera- 
tors, and premium services may be deployed if the connection 
reliability is close enough to one. Reliability is therefore a 
crucial parameter in the design and analysis of networks. 

The study of network reliability has led to a huge body of 
literature, starting with the work of Moore and Shannon [1], 
and including excellent textbooks and surveys [2], [3], [4], 
[5], [6], [7], [8]. In what follows, we consider a probabilistic 
approach, in which the network is represented by an undirected 
graph G = {V, E), where y is a set of nodes (also called 
vertices) and i? is a set of undirected edges (or links), each 
of which having a probability p„ or pe to operate correctly. 
Failures of the different constituents are assumed to occur at 
random, and to be statistically independent events. Among the 



different measures of reliability, one may single out the k- 
terminal reliability, namely the probability that a given subset 
K of k nodes {K C E) are connected. The most common 
instances are the all-terminal reliability Rel^ {K = E) and 
the two-terminal reliability Rel2(s -^ t), which deals with a 
particular connection between a source s and a destination t. 
Both of them are affine functions of each p„ and pe . 

The sheer number of possible system states, namely 
2|s|+|y|^ clearly precludes the use of an "enumeration of 
states" strategy for realistic networks, and shows that the 
final expression may be extremely cumbersome. Consequently, 
most studies have considered graphs with perfect nodes 
(Pn = 1) and edges of identical reliability p; radio broadcast 
networks have also been described by networks with perfectly 
reliable edges but imperfect nodes [9], [10]. It was shown early 
on — see for instance the discussion in [4], [11] — that the 
calculation of fc-terminal reliability is #P-hard in the general 
case, even with the restricting assumptions that (i) the graph 
is planar (ii) all nodes are perfectly reliable (iii) all edges have 
the same reliability p. All reliabilities are then expressed as a 
polynomial in p, called the reliability polynomial. 

The difficulty of the problem has stimulated many ap- 
proaches: partitioning techniques [12], sum of disjoint prod- 
ucts [13], [14], [15], [16], [17], [18], graph simplifications 
(series-parallel reductions [1], delta-wye transformations [19], 
[20], [21], factoring [22]), determination of various lower and 
upper bounds to the reliability polynomial [2], [4], [23], [24], 
[25], Monte-Carlo simulations [26], [27], [28], and ordered 
binary decision diagram (OBDD) algorithms [29], [30], [31], 
[32]. Other decomposition methods have also been proposed 
[33], [34], [35]. The reliability polynomial has been exten- 
sively studied [36], [37], [38], with the aim of finding general 
information from the structure of its coefficients [36], [37] or 
the location of its zeros in the complex plane [39]. 

In recent years, the tremendous growth of Internet traffic has 
called for a better evaluation of the reliability of connections 
in optical networks. Actual failure rates and maintenance data 
show that a proper evaluation of two-terminal reliabilities must 
put node and edge equipments on an equal footing, i.e., both 
edge (fiber links, optical amplifiers) and node (optical cross- 
connects, routers) failures must be taken into account. The 
possibility of node failure has been considered in early papers 
[9], [40]. Adaptation of algorithms to include imperfect nodes 



has been addressed [32], [41], [42], [43], [44]. In order to 
be realistic, different edge reliabilities should be used too: for 
instance, the failure rate of optical fiber links is often assumed 
to increase with their length. 

In recent works, we have shown that the two- and all- 
terminal reliabilities can be exactly calculated for recursive 
network architectures, where the underlying graphs are undi- 
rected and the edge/node reliabilities arbitrary [45], [46], [47]. 
The final expressions are products of transfer matrices, each 
element of which is a multilinear polynomial of the individual 
edge or node reliabilities constituting the 'building block' (or 
'elementary cell') of the recursive graph. 





Fig. 1. General K4 ladder. The source is always So, the destination is Sn 
otT„. 




Fig. 2. Last building block of the directed K4 ladder Edges and nodes are 
identified by their reliabilities. 

In this work, we show this general result holds for di- 
rected networks too. As an example, we calculate the two- 
terminal reliability of the directed K4 ladder displayed in 
Fig. [T] This network describes a common (nominal + backup 
paths) architecture, with additional connections between transit 
nodes enabling the so-called "local protection" policy, which 
bypasses faulty intermediate nodes and/or edges. By letting the 
individual node and edge reliabilities take arbitrary values, we 
actually do not add to the complexity of the problem but make 
the internal structure of the problem more discernible. It is then 
easier to fully exploit the recursive nature of the graph. The 
two-terminal reliability's exact expression is a product of 5 x 5 
transfer matrices, as in the undirected case [47]; consequently, 
it can also be determined for an arbitrary size (length) of the 
network. If edges have the same reliability p, and nodes the 
same reliability p, the two-terminal reliability can be expressed 
as a sum over the eigenvalues of the unique transfer matrix, 
and its generating function is a rational fraction. For large 



networks, the eigenvalue of highest modulus is, to all purposes, 
the scaling factor of the asymptotic power-law behavior. The 
determination of the failure frequency and the failure rate of 
the connection is then straightforward. Prompted by the nearly 
universal character of the Brown-Colbourn conjecture [39], 
we also address the location of complex zeros of the two- 
terminal reliability polynomial and show that they (i) may be 
quite different for directed and undirected networks (ii) exhibit 
structural transitions at various values of p. 

Our aim is (i) to give a description of the methodology 
followed in the derivation of the final results, so that re- 
searchers or engineers involved in reliability studies can use 
them, even in worksheet applications (ii) emphasize again the 
importance of algebraic structures of the underlying graphs in 
the determination of network reliability [5], [48]. 

Our paper is organized as follows. In Section [III we give 
the basic formula and methodology used in the decomposition 
method [47], which must be adapted here for directed graphs. 
In Section [Till we give the exact solution for the two-terminal 
reliability for the directed K4 ladder (with source Sq and 
destination Sn or T„). In Section |IV| we consider directed 
and undirected configurations of a special architecture. We first 
give the generating functions of the two-terminal reliabilities 
when all edges and nodes reliabilities are p and p, respectively; 
we then derive very simple analytical expressions for Rel2. 
The average failure frequency and the failure rate of the 
connection under consideration are then deduced. We then 
show that the location of the complex zeros of Rel2 differ 
in the two configurations. We conclude by proposing several 
directions in which the present results may be extended. 

II. Graph decomposition 

The purpose of our method is to simplify the graph by 
removing links of the n*^ (last) elementary cell of the network, 
namely the edges and nodes indexed by n, a procedure called 
pivotal decomposition or deletion-contraction [4]. In the case 
of undirected graphs [47], we saw that if the end terminal t 
(which can be regarded as perfect) is connected to node u 
through edge e, with respective reliabilities pu and pe, then 

Rel2(G) = {l-pe)RehiG\e)+p,p^RehiG-e) 

+Pe{l-pu)Rch{G\u), (1) 

where G \ e and G \ u are the graphs where e or u have 
been deleted, and G ■ e the graph where t and u have been 
merged through the "contraction" of e; eq. ([T]i merely sums 
probabilities of disjoint events. Here, we must adapt this 
decomposition to a directed graph. This is done very easily by 
first discarding all edges whose origin is t and then applying 
eq. ([T]i to the remaining graph. This procedure, along with 
standard series-parallel reductions, has to be repeated for the 
three secondary graphs in order to take advantage of the 
structural recursivity of the graph. After a finite number of 
such reductions, we get replicas of the original graph, albeit 
with one less elementary cell and with the (n— 1)*'' cell's edge 
and node reliabilities possibly renormalized by those of the n'^ 
cell, or set to either zero or one. In order to ensure the existence 



of a recursion relation, the graph structure must be closed 
under successive appHcations of eq. ([T]i; it may initially require 
the use of extra edges with symbolic reliabilities, so that all 
nodes of an elementary cell are connected pair-wise, even if 
such links do not exist in the graph under consideration (the 
"scaffolding principle"). At this point, a recursion hypothesis 
is needed, giving for instance Rel2(S'o -^ Sn) as a sum over 
specific polynomials in the reliabilities indexed by n; these 
are often obvious from the n — 2 value. Going from 71 — 1 
to n provides the transfer matrix linking the prefactors of 
the polynomials, because Rel2 is an affine function of each 
component reliability; the (often trivial) n = 1 case serves as 
the initial condition of the recurrence. 

III. Exact solution for the directed K4 ladder 

Let us first illustrate this method by calculating TZn = 
Rel2(S'o -^ Sn) for the directed K4 ladder shown in Fig. [Tl 
a more detailed view of the last building block is displayed 
in Fig. |2] Following the guidelines of the preceding section, 
we remove the irrelevant edges, namely those indexed by a'^, 
b'j^, and d^. The first application of eq. ([T]i is represented in 
Fig. [3] where each prefactor is put in front of the associated 
secondary graph. Note that the three secondary graphs are 
structurally identical, they differ by the actual values of the 
edge reliabilities only. Clearly, a second application of eq. ^ 
provides two-terminal reliabilities, with Sn-i or Tn-i as new 
endpoints. This should call for a similar decomposition of the 
two-terminal reliability Rel2(S'o-^T„), after which we would 
get coupled recursion relations for the two destinations 5„ and 
r„. However, this is unnecessary because the two destinations 
are identical under the permutations a„ *^ e„, c„ ^^ d„ (and 
the corresponding a'„ <-^ e^, etc.), and S'„ ^^ T„. It turns out 
that TZn may be expressed as the sum of five polynomials in 
a„, ..., Tn (see below). This is also true for Ke\2{Sa^Tn), 
which leads us to a 10 x 10 transfer matrix (the calculations 
are routinely performed by mathematical software). Because 
several lines of this matrix are identical, regrouping terms 
actually allows to limit the transfer matrix's dimension to 5, 
as in the undirected case [47]. The value of TZi, which can be 
easily calculated, leads to 



with 



7e„ = (1 0) M„ M„ 





■ Ml Ma 




(2) 



for Rcl2(5o^T„), the left vector should be (0 1 0). The 
transfer matrix Mi is given by 
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Xl = Si {ai + b'ieiTi- aib'iCiTi) , 
X2 = Si [di + b'i Ci Ti - b'i Ci di Ti) , 
2^3 = Si (fli di + ai b'i Ci Ti - Ui b'i Ci di Ti + b'i Ci Ci Ti 

-ai b'i Ci Ci Ti + (1 - a,;) b'i (1 - a) di e^ Ti) , 

(1 - ai) (1 - b'i) c'i di Ci Si Ti, 

ai (1 - b't) a (1 - dj) e'i Si Ti, 

{ci + ai bi Si - flj bi Ci Si) Ti, 



X4 
X5 

X6 

Xl 
Xs 



(ci + bi di St - bi a di Si) Tj, 



— [Ci Ci + ai bi Ci bi + ai bi di Si ai bi Ci di Si 
-ai bi a Ci Si + {1 - ai) bi (1 - q) di a Si) Ti, 
Xg = ai (1 - bi) Ci d'i (1 - Ci) Si Ti, 

a^io = a'i (1 - bi) (1 - Ci) di Ci Si Ti, 

a^ii = (1 - ai) (1 - b'i) a Si Ti, 

x\2 = (1 - 6'j) Cj (1 -di) S'jTi, 

x\-i == (1 - b'i) {aiCi- aiCtdt + CiCi- aidCi 
+ (1 - flj) (1 - Ci) di Ci) Si Ti, 

x\4 = [l-b'i) [aiCi- aiCidi + CiCi- aidCi 
+ (1 - ai) (1 - a) di Ci) Si Ti 
- {ai bi a + ai bi di - ai bi Ci di + a Ci 
-ai bi Ci Ci + {\- a.i) bi (1 - q) di a) Si Ti, 

Xi5 = {1- bi) {aiCi + aidi- aiCidi- aidCi 
+(1 - flj) (1 - Ci) di Ci) Si Ti, 

Xn = - {ci{b'^+di-b'^di) +bidi{l- Ci)) SiTi, 

xis = -{aiibi+ei-bid) + b'iei{l-ai)) SiTi, 

xi9 = (l-bi) {1-Ci) diSiTi, 

2^20 = ai (1 - bi) (1 - d) SiTi. 

For i = 0, we must set oo = do = 1 and cq ^ eg — 0. These 
formulae apply to the most general directed K4 ladder, and 
we recover the undirected case by setting a'^ — ai, b'^ — bi, 
etc. A missing edge or node is accounted for by setting the 
relevant reliability to zero, as will be seen in the following 
section. 

IV. Application 

Let us apply the results of the preceding section to the 
architecture represented in Fig.|4] The calculations are straight- 
forward, since we merely have to replace all nonexistent 
edges and nodes by zero. Does the removal of these network 
elements drastically change the previous results? Actually, no, 
even though in the directed case the dimension of the transfer 
matrix is reduced. 

A. Transfer matrices 

In both directed and undirected configurations, we must set 
Tq = Tn = and bi = (0 < i < n). This does not change the 
dimension of the transfer matrix in the undirected case, which 
remains equal to 5. However, for the directed configuration. 




b'n = (1 - bj b',.i//{c„ T„ e'„) j f b,.i//(c'„ T, ej 

dn 



+ b'„ T„ b;.i'--'bn.i 



+ b', (1-TJ b;.i-.--b„.i 




Fig. 3. First step of the decomposition. a//fe corresponds to a and 6 in parallel, and is tlierefore equal to a + b — ab. 





Fig. 4. Architecture discussed in this work ("Echelle Angele"). The source 
is So, the destination is Sn for both directed and undirected networks. 



the further simplification a^ ^ b'^ = c'l = d[ — e[ — Q is such 
that we can reduce the transfer matrix to a 3 x 3 one, namely 

CiTi CiTi CiEiT^ , (4) 

with Xi = a-i Ci (1 - di) (1 - e.^) + d^ e^ (1 - a^ - q), so that 

7^„ ^(10 0) M„ Mn^i • • • Ml Mo I I . (5) 

In many studies, edge reliabilities are considered identical 
to p in order to provide clues to the general behavior of the 
connection reliability, while nodes are viewed as perfect (i.e., 
their reliabilities are set equal to 1). In this work, we keep 
imperfect nodes with identical reliability p. We thus have two 
independent parameters to describe Rel2(S'o — > "Sn)^ allowing 
us to better distinguish the contributions of edges and nodes. 
Using only these two parameters implies that a unique transfer 
matrix needs be considered. Equations (|2]l and (|5]l show that 
the two-terminal reliability is given by the n^^ power of this 



matrix. For instance, eq. (|4|i leads to 
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P^ P 
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-p^ p^ p^ p^ [2 - Ap + p^) 



(6) 



B. Generating functions 



Because of the intrinsic recursion relation between succes- 
sive powers of M{p, p), a similar one should hold for the two- 
terminal reliability. The generating function formalism [49] 
is a useful way to store all necessary information in a very 
concise manner It is defined by 



Q{^) 



n=0 



Rel2(5o^5„)z", 



(7) 



Its calculation is straightforward. Here, G{z) is necessarily a 
rational fraction of z, namely 



Giz) 



Af{z)_ 
V{z)- 



(8) 



Using eqs. Q and (|6]l (for the directed configuration), we 
compute the first values of Rel2(S'o — > Sn) for, say, n equal 
to 2, 3, 4, etc., in order to obtain the first terms of the expansion 
in z of Q{z). We then multiply this truncated expression by 
the characteristic polynomial of the transfer matrix taken at 
1/z. The series expansion of this product in the vicinity of 
z = leads to M{z). No determination of eigenvalues or 
eigenvectors of the transfer matrix is needed. 

Having described the method, we limit ourselves to the final, 
simplified expressions, which should be used for n > 2, even 
though they are valid for n ~ \ (both give the first-order term 
p p^ z, corresponding to the two-terminal reliability between 
the two nodes 5o and Si, connected by a single edge). 



1) undirected case: 
P 



Mu{z) = 



[l~p (2 + 2p-6p2 + 3p3) p2 



2{l-p) 

+2 (1 - pf (2 - p) p^ {1-P + p") p^ z^) , (9) 

Vuiz) = 1-pp {2 + 2pp-6p^p + 3p^p) z 

+2 (l-p)p3p3 (2-3p-2p/? + 6pV 

-Ap^p + p'^p) z^ 



5 ^3 



-4(l-p)^ i2~p)p' (l-p) p^'z' 

2) directed case: 

Afd{z) = ^ {l-p^p^2-ip + p^)z), 

Vd{z) = 1-pp {2 + 2pp-4p'^p + p^p) z 
+2p3p3 (i_p) {2-p)z\ 



(10) 



(11) 



(12) 



C. Analytical expressions of the two-terminal reliabilities 

The two-terminal reliabilities are derived from the partial 
fraction decomposition of the associated generating functions 
Gu{z) and Gd{z), because the eigenvalues of the transfer 
matrix are the inverses of the generating function poles. 

1) undirected case: Equations (l9moli can be further sim- 
plified for perfect nodes {p — 1), because T>u{z) is then of 
degree 2 in z, from which we find 



nl^Hp^ 


I;n>2)^a+(C1"') + a_ (ci"^ 
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An) 


1 (2 + 2p- 6p2 + 3p^ ± VA^ , 




(14) 
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B±{l-Ap + 2p^) VA 
iil-pfil-p + p^fVA' 
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A - 


4 - 8p + 36 p2 - 100/ + 128/ 








-76/ + 17/ , 




(16) 


B = 


2 - lOp + 38 p2 - 59/ + 40/ - 


10 p5 


.(17) 



For p < 1, we would get a sum over three eigenvalues, the an- 
alytical expression of which would only be more cumbersome 
(it will not be given here). 

2) directed case: Even if p ^ 1, Vdiz) is of degree 2 in 
z. We have again two eigenvalues, so that 



7^W(n>2)=a+ C 



.(d) 



a- C 



.(d) 
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Ad) 



a± 
A' 



^ (2+pp{2-ip + p^)±VA' 
^^2-pp(2-4p + / 



4 V VA' 

4-4pp(2-2p + /) 
+/p^(2-4p + p2)2. 



(18) 

(19) 
(20) 

(21) 



(^ and Cj_ are displayed in Fig. \5\ for perfect nodes. 
Obviously, C+ is nearly always much larger than (-, especially 
when p lies in the vicinity of 1 . (^V" and (^\ are nearly equal 



over the whole range < p < 1, while the second eigenvalue 
is much larger in the directed case. Still, as n grows, the 
contribution of the second eigenvalue should vanish so that 
the two-terminal reliability exhibits an asymptotic power-law 
behavior TZn oc C", the scaling factor being C+, the eigenvalue 
of largest modulus. Even for n w 10, this asymptotic limit 
would already be a good approximation. 




Fig. 5. Variation of f+ (full line) and (— (dashed line) with p for the directed 
(red curve) and undirected (purple curve) architectures of Fig. |4] with perfect 
nodes. 



D. Average failure frequency and failure rate 

Steady-state system availability A and failure frequency j7 
are important performance measures of a repairable system, 
from which other key parameters such as the mean time 
between failures, average failure rate, Birnbaum importance, 
etc. may be deduced [6], [7], [8]. The first calculations of 
steady-state failure frequencies were based on the inclusion- 
exclusion principle, with adequate failure and repair rates 
(more generally, the inverses of the mean down and up times) 
were attributed to each term of the expansion [50], [51]. 
Several papers have since provided a few simple recipes, 
describing how F and the system failure rate A == V/A can then 
be derived [52], [53], [54], [55]. All these formal calculations 
boil down to a simple fact: the failure frequency may be 
derived from A by the application of a linear differential 
operator [56], [57], [58]: 



t^ = y^ AiP, 



dA 
dpi 



Z^ 



fJ-iQi 



dU_ 
dqi 



(22) 



where [/ = I — A is the total unavailability, pi the availability, 
Qi = 1 — Pi the unavailability, A^ the failure rate, and pi 
the repair rate of equipment i. The expressions obtained in 
the preceding sections {A = TZn) make such calculations 
straightforward by the application of the linear differential 
operator to each transfer matrix. Good estimates of what 
happens for large networks (n ^ 1) may be obtained by 
considering that all links have availability p and failure rate A 
(assuming nodes are perfect, to keep the discussion simple). 
If TZn ~ fl+ C+' the average failure rate A„ is then given by 



^n — Z^ — 
7? 



A 



d 111 a+ 
d\np 



glnC+ 
dlnp 



(23) 




Fig. 6. Variation of dln(^^/dlnp with p for the directed (red curve) and 
undirected (purple curve) architecture of Fig. |4] with perfect nodes. 



The "effective" network can therefore be seen as a series 
network of n components, having each an availability equal to 



C+ and a failure rate 



d\np 



X. From Fig.|5]we could expect that 



the asymptotic failure rate for our case study would be roughly 
similar for directed and undirected configurations, over the 
whole range < p < 1. This is indeed observed in Fig. |6l 

E. Zeros of the reliability polynomials 

The structure of the different reliability polynomials may 
be understood by studying the locations of their zeros in the 
complex plane. Such a study has been fruitfully performed for 
chromatic polynomials [59], [60], [61]. In reliability studies, 
a search for general properties of the all-terminal reliability 
polynomial RelA(p) [36], [37], [38] has brought the Brown- 
Colbourn conjecture [39], according to which all the zeros 
could be found in the region |1 — p| < 1. Although valid 
for series-parallel graphs, this remarkable conjecture does not 
strictly hold in the general case (but not by far) [62] . Rel^ {p) 
is linked to the Tutte polynomial, a graph invariant. It has 
also been studied extensively by Chang and Shrock for various 
recursive families of graphs [63], who give the limiting curves 
where all the zeros converge. 

As n grows, the number of complex zeros of the reliability 
polynomial increases. Because of the matrix transfer property, 
there is a recursion relation between reliability polynomials 
corresponding to successive values of n. The general problem 
has been treated by Beraha, Kahane, and Weiss [64]. It may be 
understood in the following, simplifying way: if the reliability 
polynomial is of the form X]j "^j '^j(p)" (where Xi are the 
eigenvalues of the recurrence), then at large n, only the two 
eigenvalues of greater modules, say Ai and A2, will prevail, 
so that the reliability polynomial will vanish when |Ai(p)| = 
|A2(p)|. This equality defines a set of curves in the complex 
plane, where all zeros should accumulate in the n ^r 00 limit. 
A detailed discussion of the convergence to the limiting curves 
may be found in [61]. 

The location of the zeros of Rel2 [p) in the complex plane is 
also worth investigating, even though Rcl2(p) is not a graph 
invariant. The new twist lies in the extra parameter at our 
disposal, the node reliability p, which has a deep impact on 



the curves to which the zeros of Rel2(p) converge as n ^ 00. 
Actually, structural changes occur at critical values of p, which 
can be deduced from the expressions of the eigenvalues. 

We have displayed in Figs. |7] and |8] the location of com- 
plex zeros of the two-terminal reliability polynomials from 
our case study (the numerical values have been obtained 
by using Mathematica). The structures of these zeros are 
quite different for the directed and undirected cases. In both 
configurations, parts of the real positive axis are indeed the 
limiting "curve". It turns out that in the undirected case, there 
is no such real line segment for exactly pc^ = | (the sole 
intersection of the right-most curve with the real axis lies at 
p — |). We shall not detail the other critical values for p. 
In the directed case, the segment gets closer and "punctures" 
the curve on the right half-plane for pca ~ 0.51242 (a root 
of a polynomial of degree 10). As p further decreases, both 
structures look like circles plus an extra line segment, which 
expands as p^^/"^. We limit ourselves to the asymptotic values 
of IpcIicIcI (the circle) and p± (the endpoints of the segment 
on the real positive axis). Note the expansion rates are distinct, 
too. 

1} undirected case: 
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2) directed case: 
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V. Conclusion and outlook 

The two-terminal reliability of directed networks may also 
be expressed by a product of transfer matrices, in which each 
edge and node reliability is exactly taken into account. The 
size of the transfer matrix should of course increase with the 
network's "width". This general result could be extended to 
the all-terminal reliability with nonuniform links [65]. We can 
now go beyond series-parallel simplifications and look for new 
(wider) families of exactly solvable, meshed architectures that 
may be useful for general reliability studies (as building blocks 
for more complex networks), for the enumeration of self- 
avoiding walks on lattices, and for directed percolation with 
imperfect bonds and sites. Since the true generating function 
is itself a rational fraction, Pade approximants could provide 
efficient upper or lower bounds for these studies. Moreover, 
individual reliabilities can be viewed as average values of 
random variables. Having access to each edge or node allows 
the introduction of disorder or correlations in calculations [66]. 
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Fig. 7. Location of complex zeros for the undirected ladder, with n = 100 
and p equal to 1, 0.9, and 0.1 (from top to bottom). 



Fig. 8. Location of complex zeros for the directed ladder, with n = 50 and 
p equal to 1. 0.9, and 0.1 (from top to bottom). 
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